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We investigate the interplay between long range and direct photo-assisted transport in a triple 
quantum dot chain where local ac voltages are applied to the outer dots. We propose the phase 
difference between the two ac voltages as external parameter, which can be easily tuned to manipu¬ 
late the current characteristics. For gate voltages in phase opposition we find quantum destructive 
interferences analogous to the interferences in closed loop undriven triple dots. As the voltages 
oscillate in phase, interferences between multiple paths give rise to dark states. Those totally cancel 
the current, and could be experimentally resolved. 

PACS numbers: 73.63.Kv, 75.10.Jm, 85.35.Be, 85.35.Gv 


I. INTRODUCTION 

A system that is driven non-adiabatically through the 
avoided crossing of two states undergoes a transition^^^. 
The probability of the transition depends on the param¬ 
eters of the driving and the splitting at the crossing. 
The latest is given by the coupling between the diabatic 
states. Repeating the passing through the crossing in¬ 
troduces different paths to end in a given state, which 
gives rise to constructive interference. The control of 
this mechanism in solid state qubits has become a stan¬ 
dard tool in the manipulation of quantum states^^^, the 
generation of entanglement®, or the measurement of the 
qubit coherence timescales®. 

In periodically driven quantum dot systems, this effect 
is measured as photon-assisted tunneling resonances^®. 
An electron is hence delocalized between tunnel coupled 
quantum dots when the detuning of their energy levels 
is a multiple of the driving frequency nhuj^^. The tun¬ 
nel coupling is renormalized by the ac field by a factor 
which depends on the amplitude and frequency of the 
driving^®’^®. Recently, striking electron spin resonance 
measurements in quantum dot systems®^ have been in¬ 
terpreted in terms of multilevel crossings®®. Three-level 
crossings may also lead to peculiar phenomena such as 
dark resonances®®’®^. 

Triple quantum dots (TQDs) are ideal systems for in¬ 
vestigation of such processes. On one hand, the spacial 
separation of three states®®’®®, one in each dot (L, C and 
R) makes it possible to manipulate them individually by 
means of gate voltages®®. Hence, different drivings can be 
applied to the different levels by applying localized time- 
dependent gate voltages to each quantum dot®®. Thus, 
not only the amplitude and frequency of the driving®®’®®, 
but also phase differences®® become important. 

On the other hand, the tunnel coupling between all 
three states can be tuned, also between those that are 
not directly coupled. Indeed, long-range transport be- 
tweeen the edge dots of a linear TQD has been very re¬ 
cently detected®^^®®. During these higher-order (cotun¬ 
neling) transitions, the center dot is only virtually occu¬ 


pied. Hence they involve the direct transfer of a charge or 
a spin qubit betweeen distant sites, avoiding decoherence 
and relaxation in the intermediate region. 

By applying sinusoidal signals to the outer dot gates, 
the system can be driven through anticrossings among 
the three states, performing Landau-Zener transitions: 
L-C, G-R and L-R, cf. Fig. 1(b). Note that the long- 
range L-R transition is paralell to the direct-tunneling L- 



FIG. 1: (a) TQD in series connected to leads. Two ac volt¬ 

ages are applied to the outer dots, with a phase difference 
4>. (b) Time evolution of the energy levels for the left, center 
and right dots for </> = tt, showing the different crossings me¬ 
diated by direct (L-C and C-R) and virtual tunneling (L-R). 
(c) For (f) = Ti the driving induces resonant transitions (ei¬ 
ther direct, or virtual, T(()^") between all the levels. The 
system can thus be mapped to an undriven triangular TQD 
depicted below, (d) The case where the edge dot levels os¬ 
cillate in phase {(j> = 0) is equivalent to that where only the 
center dot is driven. Transport is then governed by sidebands 
in the center dot. 
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C, C-R trajectory. The system then behaves as a combi¬ 
nation of coupled interferometers. The interference pat¬ 
terns, coming from real or virtual paths, can be detected 
by coupling the system to source and drain contacts and 
measuring the current that passes through it along a volt¬ 
age bias eV, as sketched in Fig. 1. The versatility of this 
system makes it suitable for investigating many differ¬ 
ent configurations. Here we will consider the case where 
the two drivings are equal in amplitude and frequency, 
ki(i) = (14c/2) cos(a;t -I- (j)i), but they have a phase dif¬ 
ference 0 = (/>R — ^L- As we show below, the dependence 
on the phase is of paramount importance and it has not 
yet been addressed. 

We concentrate on two cases, when the dots oscillate 
with opposite ((/> = tt) or the same phase {(j) = 0). In 
the first case, </) = tt, crossings of all three states occur. 
Direct tunneling transitions between left and center and 
center and right take place when the corresponding en¬ 
ergy levels satisfy £c — El ~ mhui and £r — £c ~ m'huj. 
Importantly, even if they are not directly coupled, the 
crossing of the outer dot levels induces a long-range tun¬ 
neling transition. It is mediated by the virtual occupa¬ 
tion of the center dot, when er — en ~ nhuj, regardless 
of the energy Thus, different paths are possible 

for left to right transport [see Fig. 1(c)] remarkably lead¬ 
ing to the interference of real and virtual (cotunneling) 
transitions. 

Differently if the two dots oscillate in phase {(j) = 0), 
only L-C and C-R crossings are possible. The system 
is therefore equivalent to the one where only the cen¬ 
ter dot is driven. Hence £c develops sidebands sepa¬ 
rated by an energy hoj, as shown in Fig. 1(d). Resonant 
transitions involving the sidebands is the main transport 
mechanism. However, long-range tunneling is also possi¬ 
ble if the outer levels have the same energy. The virtual 
occupation of sidebands with either positive or negative 
detuning causes unexpected interferences, similar to two 
path interferometers®®’®^. 

The paper is organized as follows: we present the 
model and numerical results in Sec. H. An analytical ap¬ 
proach is performed in Sec. HI which allows us to give 
an interpretation of the main features in Sec. IV. A sum¬ 
mary is given in Sec. V. 


II. MODEL 


We assume the system to be in the Coulomb block¬ 
ade regime, where only one electron is allowed in the 
TQD at a time. In such a case, the spin of the elec¬ 
tron does not play any role and we can ignore it. The 
states of our Hamiltonian are defined in the basis {\L), 
jC), |i?), jo)}, where |i) represents one electron in the 
i-dot, and | 0 ) corresponds to the case where the system 
is empty. Furthermore, we consider a large bias voltage 
such that electrons flow unidirectionally from the left to 
the right reservoir. We write the Hamiltonian of the sys¬ 


tem H{t) = Hq + iFiead + Hint + Hnc{t), where 

Hq = He+Hr = ^EicJci-b^(n,j+i4c*+l -f H.C.) 

i i 

( 1 ) 


describes the TQD with energy levels Si localized in each 
dot and (nearest neighbor) interdot tunnel couplings Tij. 
The term Hiead = 'Yliik ^ikdli.dik represents the leads, and 
Hint = + H-C. is the system-lead interac¬ 

tion Hamiltonian, given by the coupling 7 ;. The driving 
Hnc(t) = ^cts on the edge dots /=L,R. 

In order to treat the explicit time dependence, it is 
convenient to perform a unitary transformation H'(t) = 


H{t) — iHdt U{t), where 


U{t) = exp 


dtiHncih 


(2) 


It leaves the diagonal terms time independent. The sys¬ 
tem Hamiltonian transforms as Hq^I) = H^ + H'^.^ with 


00 

H'r = ^ -I- t'^^e -I- H.c. (3) 

U— — 00 


where the hopping terms = tlcJ v {Vac/ 

= 'TcRJv{Vnc/2huj)e~^''^ are renormalized by the 
i/-th order Bessel functions Juix). In what follows, we 
will simply write J,y, unless its argument is different from 
Vnc/(2hui). In the large bias regime we are considering, 
eV ^ hu!, Vac) tbe effect of the transformation on the 
tunneling term Hint can be disregarded due to the nor¬ 
malization of Bessel functions. 

In the weak coupling limit and assuming the Born- 
Markov secular approximation®®, the quantum master 
equation for the reduced density matrix of the TQD be- 

SQ 

comes : 

m = '-[H,+H{{t),p{t)] + iCr-U)pit) (4) 

where 

(m1^) “ ^ ^ (r nk ^kn^Pmn{i')^7nn') (5) 

k^n 

“ 2 I ^ ^ ^ ^ ^km j Pmn(l 

\k^n k^m j 


In the large bias regime, only two rates participate, 
Tlo = Tl and Fqr = Fr, with F; = {2tt/ h)'Di\^i\^, and 
Vi being the density of states in the lead to which dot I is 
coupled. The steady state occupations are obtained 
by numerically integrating in time Eq. (4). With those, 
we get the stationary current, given by / = gErpIJj^. 






3 


I/eTj, 





■(d) J 


i 



/HO 6 6.4 6.S 

v^/(2nuj) 

0.3r 


-7 -6.5 -6 

ec/?i“ 



0.5 1 1.5 

yac/(2?iw) 


FIG. 2: (a) Current through a TQD by tuning ec and Vac. 
ER = EL = 0, w = S.Stlc, tlc = TCR and ij) = it. Superim¬ 
posed to the well-known Landau-Zener-Stiickelberg pattern, 
we find exact cancellations of the current that are attributed 
to the contribution of virtual tunneling transitions. The green 
dashed lines correspond to the zeros of Jo{Va,c/huj), cf. Eq. 
(10). (b) and (c) show the current for two cuts along the 

resonance ec = —Shuj [marked as red lines in (a)]. While 
in (c) the current vanishes coinciding with JoiVaA^/Tiw) = 0, 
a double suppression appears in (b) around the first zero of 
JiiVa.c/'ihuj)- The complicated double-minimum interference 
patterns are of the same nature as the ones observed in an 
undriven triangular configuration. By means of the RWA one 
can derive an analytical condition Eq. (8) [blue dashed lines 
in the inset of (a)] for the dark state Eq. (9). (d) current 
for I4c/(2/ia;) = 3.51 in the cotunnel regime, where Fano-like 
resonances are observed. 


A. Numerical observations 


The current for the different configurations, (jj = tt and 
(/) = 0, is plotted in Figs. 2 and 3 respectively, as func¬ 
tions of the driving amplitude, t4c, and the detuning of 
the centre dot level ec from ep = cr = 0. In both con¬ 
figurations, we can clearly distinguish two regions by the 
ratio a = V"ac/(2|ec|), which marks the onset of photon- 
assisted transitions. The current shows very different fea¬ 
tures in each region, which we detail below. 

For (j) = TT, the current is cancelled for certain driv¬ 
ing amplitudes in the regime a < 1 where the center dot 
is not excited, see horizontal lines in Fig. 2. For larger 
amplitudes, we observe the expected photon-assisted res¬ 
onances around ec = rnhoj which are suppressed at the 
f-th zeros Zi^m of More surprisingly, we find addi¬ 



ec/Tvzi 


FIG. 3: Gurrent through the TQD as a function of the de¬ 
tuning ec and 14 c for <)> = 0, with er = el = 0, ta = 3.3rLc, 
Tlc = tcr. Narrow antiresonances are due to destructive in¬ 
terference between sidebands with positive and negative de¬ 
tuning with respect to the L-R resonance. The dashed lines 
within the white regions correspond to the fulfillment of condi¬ 
tion (12) for dark states, see text below. Right panel: current 
at fixed ec = —3hu) (red vertical line). The dips correspond to 
the zeros of J3. Bottom panel: current cut along 14c = Shiz 
(horizontal blue line) showing a symmetric configuration of 
zeros which comes from the interference between sidebands. 


tional off-resonance pairs of suppressions of the current 
in between zeros of consecutive resonances, cf. inset in 
Fig. 2. Remarkably, the cancellation of the current is 
exact. These sharp features are a manifestation of the 
interference of real and virtual transitions, as we discuss 
in the next section. 

Very differently, for ^ = 0, no significant features are 
observed below 0 = 1, cf. Fig. 3. For larger amplitudes, 
we find narrow off-resonance current suppressions all the 
way between the f-th zeros Zi^m of increasing/decreasing 
order, e.g. current is cancelled between Zi^m and 
cf. Fig. 3. These features come from the destructive in¬ 
terference between sidebands with positive and negative 
detuning with respect to the L-R resonance. 


III. ANALYTICAL APPROACH 

We perform an analytical treatment consisting on a 
perturbative expansion of the evolution operator in the 
inter-dot tunneling, and a rotating wave approximation 
(RWA) close to the relevant resonances. Direct tunneling 
between nearby dots (L-C and C-R) is contained in the 
first order of the expansion. The long-range tunneling 
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between the left and right dots is only captured in the 
second order. 

We are interested in the configurations where both 
real and virtual transitions coexist. This occurs close 
to the double resonance conditions ec — £l ~ mhuj and 
£r — £l ~ nhoj, as depicted in Fig. 1. Close to these, 
one can transform the Hamiltonian to the rotating frame 
and neglect the fast oscillating terms of the Hamiltonian. 
With this RWA we obtain a time independent Hamilto¬ 
nian: 




i—m 

^LC 

T™/(</.) \ 


i—m 

^CL 

~m+n 1 A m,n 

tcR-^W 


V 

iRC-”(</') 

4+Arr” / 


where ef = £i + nhuj. Xlrr = r^^Xm+n 

and = — (A™l" -|- A^^ ) are the level shifts coming 

from the second order tunnel contribution, with Xp = 
— Jp)Jvl{(p— v)fujj\. More importantly. 


TlTif W = t'lctcr ^ Ju 


^Ju+n Jm-\-2ri)£ 

2 (z/ — m — n)Huj 

( j piv(j> _ T 

^m—n^ J ^ 

2{v — m)hio 


( 7 ) 


describes the second-order photon-assisted tunneling be¬ 
tween the outermost dots. 


IV. INTERPRETATION 


We emphasize the role of the phase difference in the 
off-diagonal elements in Eq. (6). As we discuss in the 
following, the origin of the destructive interference pat¬ 
terns is very different in the two configurations of interest 
here, (/) = 0 , tt. 


A. Opposite phase drivings, (j) = n 


where: . Eq. (8) is the exact condi¬ 

tion for the existence of a dark-state: 

ITds) = N \L) - |C)) , (9) 

where N is a normalization constant. Note that it con¬ 
tains no contribution of the state |R), which is coupled 
to the drain lead. Therefore, no current flows through 
the system. 

The dark state condition (8) is fulfilled for all the dips 
observed in Fig. 2, showing that the linear driven system 
behaves as an undriven triangular one around every z, „j- 
The condition (8) is represented by dashed lines in the 
inset of Fig. 2 with a remarkable agreement. We can 
thus interpret these dips as the destructive interference 
between first- and higher-order tunneling transitions. 

When the ac amplitude is small, a < 1, we have 
Jm "C 1, for m > 1. Thus, in this regime the effec¬ 
tive direct couplings ^CLi^Cr" negligible. This is not 
the case for the long-range coupling. Using the addition 
theorem for Bessel functions, it can be approximated to 
^lr” ~'''i^cTCRJn{^) /mhuj. Note that the argument 
of the Bessel function is doubled and therefore it is not 
necessarily negligible. Hence, transport is dominated by 
cotunneling and Eq. (6) can be mapped to an effective 
two level system^°: 


H, 


RWA 


(tt) 


(^L,n 



TLCT'CR ^ 

mhio 


( 10 ) 


with the Pauli matrices (t^. With this model, we can 
interpret the different cancellations of the current marked 
by horizontal dashed lines in Fig. 2. They are due to 
the zeros of Jo{^), therefore representing a long-range 
analog of coherent destruction of tunneling^^. 

In this regime {a < 1), the center dot remains un¬ 
coupled from the background cotunneling transport. By 
increasing t4c, the resonant photon-assisted L-C and C- 
R transitions are activated. The center dot can then 
be considered as a localized state coupled to a contin¬ 
uum leading to asymmetric Fano-like resonances^^’^®, see 
Fig. 2(d). 


Let us start by the case (j) = tt, plotted in Fig. 2. The 
double dip detailed in the inset can be understood di¬ 
rectly by looking at the structure of the effective Hamil¬ 
tonian in Eq. (6). The long-range coupling of left and 
right dots makes it analogous with a TQD in a closed- 
loop triangular configuration. Such a system is known to 
lead to dark states when different paths around the trian¬ 
gle interfere destructively^^. This occurs when particular 
relations between couplings and detuning are met^^. In 
our case, with cl.o = cr.o = ec,m, we find the following 
condition for a dark state: 

(<(.) = tCc (Tlr°(</>)^ - tcm") ( 8 ) 


B. In-phase drivings, (j> = 0 

We now consider the configuration with (j) = 0. In 
this case Eq. (8) is only fulfilled for the trivial solution 
when all the couplings are zero. As the driving does 
not induce the left and right levels to cross, cotunneling 
is only effective when el = £r, i-e., n = 0. In that 
case, the two levels cross the center one simultaneously. 
This situation can be mapped into having the driving 
in the center dot and sidebands at energies ec + i^huj, cf. 
Fig. 1(d). If the center dot is detuned ec = mfiw + A, the 
different sidebands will act as parallel channels for the 
cotunneling transition. The effective Hamiltonian then 
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reads: 

OO 

~ +'^lct’cr X! (11) 

y——oo 

with the sideband-dependent couplings = 

~ fn)huj — A], Note that they can be tuned by 
the driving amplitude and the detuning. If the condition 

ra OO 

^ r!,,^(A) = - ^ 11.,^(A) (12) 

1/— — 00 


is met, sidebands with positive {v > m) and negative 
(y < m) detuning will destructively interfere. This leads 
to the narrow cancellations of the current shown in Fig. 3. 
The condition (12) is represented there by dashed lines. 
We note that the agreement with the full numerics is ex¬ 
cellent in all the configuration space. The system then 
behaves as an ac driven interferometer of multiple tran¬ 
sitions. 


C. Arbitrary phase difference 

The effective Hamiltonian ( 6 ) as well as the dark-state 
condition (8) are in general valid for any value of 4>- In 
this paper we focus on the simplest cases 4> = { 0 , tt} 
where all the couplings are real. For an arbitrary phase 
difference, and are complex in general. 

Therefore, Eq. ( 8 ) has to be fulfilled for both the real and 
the imaginary parts. This two conditions make the ana¬ 
lytical study of these phases much more involved, which 
is out of the purpose of this paper. We however present 
numerical results of the phase-difference dependence in 
Fig. 4 for two configurations: in Fig. 4(a) the center dot 
is resonant, ec ~ El.r = Shut, whereas in Fig. 4(b) it is 
slightly off resonance. 

As mentioned above, in resonance conditions exact 
cancellations of the current appear for (j) = rtt, n € Z 
if Eq. ( 8 ) or Eq. (12) is met for n odd or even, respec¬ 
tively. Nevertheless, we observe in Fig. 4(a) a substan¬ 
tial current drop at Vac/(21iw) « 6.3 and 9.8 for any 0, 
see Fig. 4(a). This feature is related to the zeros of the 
corresponding Bessel function (in this case the first and 
second zeros of J 3 , 21,3 and ^ 2 , 3 ). In the off-resonant con¬ 
figuration the current drop at those 14 c is not observed, 
since more sidebands are taking part in transport, cf. 
Fig. 4(b). Thus the vanishing current at this conditions 
for (j) = riTT can only be attributed to interference effects. 
In the previous sections IV A and IV B, we have discussed 
that the dark states a,t (p = 0 and 0 = tt are due to dif¬ 
ferent mechanisms. This becomes clear as they behave 
differently with a change of p. Those at </> = 0 are much 
narrower because the central symmetry depicted in Fig. 
1 (d) is broken for ^ ^ 0 . 

Other dark states can appear for specific (p mr) 
which now depend in all the other parameters. This is 
the case of features observed around 14 c/( 2 Iia;) « 2.1 



FIG. 4: Current through the TQD numerically calculated 
with Eq. (4) by tuning p and Vac. In both plots £r = er = 0, 
LU = S.Stlc and tlc = "tcr. (a) The central dot is in resonance 
with the third side-band of the two outer dots, ec = el — Sfica. 
A cut at 0 = TT corresponds to Fig. 2(b) and (c), and at 
p = 0 to the right vertical cut of Fig. 3. (b) Same as in (a) 
but with the central dot out of resonance with the outer dots, 
ec = Cl — 3.2huj. 


for p « ± 0 . 687 r in Fig. 4, which fulfill the dark-state 
condition ( 8 ). 

For smaller amplitudes, the system is in the cotun¬ 
neling regime. We find that the coherent destruction of 
cotunneling —which at 0 = tt occurs for the condition 
JoiVas/^) = 0, cf Sec. IV A— is robust and survives in 
a wide range of phase differences. 


V. CONCLUSIONS 

In summary, we predict quantum interferences that 
depend in a non-trivial way on the phase difference of 
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the locally applied drivings. For gate voltages in phase 
opposition, we find destructive interferences between 
direct and long-range transitions which are analogous to 
dark states in close loop undriven triple dot molecules. 
As the edge dots levels oscillate in phase, quantum 
paths mediated by positive and negative detuned side¬ 
bands interfere leading to multiple dark states in the 
Landau-Zener-Stiickelberg pattern. These destructive 
interferences can be experimentally detected as they 
are of the same nature as long-range current resonances 
which have been unambiguously observed. We propose 
a transport configuration, where all parameters are 
experimentally controllable, in which these features can 
be measured as cancellations of the current. This is 


particularly accessible in quantum dot architectures 
which are within experimental reach^'^’^®’^® for both 
electric or magnetic field drivings. 
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